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ABSTRACT 

Some progenitor models of gamma-ray bursts (GRBs) (e.g., collapsars) may produce 
anisotropic jets in which the energy per unit solid angle is a power-law function of the 
angle (oc 6~~ k ). We calculate light curves and spectra for GRB afterglows when such 
jets expand either in the interstellar medium or in the wind medium. In particular, 
we take into account two kinds of wind: one (n oc r~ 3 / 2 ) possibly from a typical red 
supergiant star and another (n oc r~ 2 ) possibly from a Wolf-Rayet star. We find that 
in each type of medium, one break appears in the late-time afterglow light curve for 
small k but becomes weaker and smoother as k increases. When k > 2, the break seems 
to disappear but the afterglow decays rapidly. Thus, one expects that the emission 
from expanding, highly anisotropic jets provides a plausible explanation for some 
rapidly fading afteglows whose light curves have no break. We also present good fits 
to the optical afterglow light curve of GRB 991208. Finally, we argue that this burst 
might arise from a highly anisotropic jet expanding in the wind (n oc r -3 / 2 ) from a red 
supergiant to interpret the observed radio-to-optical-band afterglow data (spectrum 
and light curve). 
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Introduction 



Although the study of gamma-ray bursts (GRBs) has been revolutionized due to observations 
of their multiwavelength afterglows and particularly due to determinations of their redshifts in the 
past few years, the question whether the GRB emission is spherical or jet-like has been unsolved. 
This question has important implications on almost all aspects of the GRB phenomenon, e.g., the 
total energy that is released in an explosion, the event rate, the physical ejection mechanism and 
the afterglow decay rate (for a review see Sari 2000). 

GRB 990123, the most energetic GRB well-studied, has an isotropic gamma-ray release of 
~ 3.4 x 10 54 ergs, which corresponds to the rest-mass energy of ~ 1.9M (Kulkarni et al. 1999). 
From such an energetics, one deduced for the first time that the GRB emission may be highly 
collimated in this case with a half opening angle of 9 m < 0.2 so that the intrinsic explosive energy 
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could be reduced to ~ 10 51 - 10 52 ergs, which is still consistent with the stellar death models, e.g., 
collapsars (MacFadyen, Woosley & Heger 2000). Another well-known feature of this burst is that 
its R-band afterglow light curve began to steepen at about day 2 after the burst (Kulkarni et al. 
1999; Castro-Tirado et al. 1999; Fruchter et al. 1999). This steepening was first argued to be due 
to the possibility that a jet has evolved from the spherical-like phase to the sideways-expansion 
phase (Rhoads 1999; Sari, Piran & Halpern 1999) or that we have observed the edge of the jet 
(Panaitescu & Meszaros 1999; Meszaros & Rees 1999). Alternative explanation for this steepening 
was subsequently proposed to be that a shock expanding in a dense medium has undergone the 
transition from the relativistic phase to the non-relativistic phase (Dai & Lu 1999, 2000). Since 
then, the dynamics of expanding jets has been numerically studied by many authors (Panaitescu 
& Meszaros 1999; Moderski, Sikora k, Bulik 2000; Huang et al. 2000a, b; Kumar & Panaitescu 
2000; Gou et al. 2000; Wei & Lu 2000). Observationally, one marked break has also been observed 
to appear in the optical afterglow light curves of GRB 990510 (Harrison et al. 1999; Stanek et 
al. 1999), GRB 991216 (Halpern et al. 2000; Sagar et al. 2000a) and GRB 000301C (Rhoads & 
Fruchter 2000; Masetti et al. 2000; Jensen et al. 2000; Sagar et al. 2000b; Berger et al. 2000). 
Holland et al. (2000) recently collected, re-analyzed and explained all of the published photometry 
for GRB 990123 and GRB 990510. All these studies have assumed that both the energy per unit 
solid angle and the bulk Lorentz factor at any angle within jets are independent of the angle. 
Hereafter we refer to such jets as isotropic ones. 

In fact, both the energy per unit solid angle (dE/dfl) and the bulk Lorentz factor within a 
realistic jet may strongly depend on the angle 6. For example, dE/dQ oc 6~ k where k ~ 3 within 
the energetic jets expected in the collapsar model for GRBs of MacFadyen et al. (2000). We call 
such jets anisotropic ones. Salmonson (2000) argued that if GRBs arise from anisotropic jets, then 
the range of the viewing angle (viz., the angle between the line of sight and the jet axis) produces 
a range in the observed properties of GRBs, i.e. the lag-luminosity relationship. In this paper, 
we calculate the emission from anisotropic jets expanding both in the interstellar medium (ISM) 
and in the wind medium. In particular, we consider two kinds of wind: one (n oc r~ 3 / 2 ) possibly 
from a typical red supergiant star and another (n oc r~ 2 ) possibly from a Wolf-Rayet star. To our 
knowledge, this is the first numerical work about evolution of anisotropic jets and their emission. 
A preliminary analysis of afterglows from anisotropic jets was made by Meszaros, Rees & Wijers 
(1998), but our results are different from their analytical ones because we take into account 
the light travel effects related to different sub-jets within the jet. We also carry out a detailed 
modelling of the optical afterglow flux data of GRB 991208 given by Sagar et al. (2000a) and 
Castro-Tirado et al. (2000). Finally, we argue that this burst might arise from a highly anisotropic 
jet expanding in the wind from a red supergiant to explain the observed radio-to-optical-band 
afterglow data. 
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2. The Model 

An energetic anisotropic jet is supposed to be ejected in an explosion. We assume that such 
a jet is adiabatic, i.e., the radiative energy is a negligible fraction of the total energy of the jet. 
This assumption is valid if the energy density of the electrons accelerated by a shock, produced by 
the interaction of the jet with its surrounding medium, is a small fraction e e of the total energy 
density of the shocked medium or if most of the electrons are adiabatic, i.e., their radiative cooling 
timescale is larger than that of the jet expansion (Dai, Huang & Lu 1999; Bottcher & Dermer 
2000). 



2.1. Dynamics 

The central question in the next study is how the Lorentz factor 7 at angle 9 within the jet 
evolves with the observer time t because all other quantities that appear in the observed flux are 
functions of 7 and of the jet radius and medium density. To obtain 7, we assume that the kinetic 
energy per unit solid angle within the jet has the following distribution: 

m ._dE j so H0<6, 

where 9 m is the half opening angle of the jet. The ejected mass per unit solid angle, M e j(0) = M , 
is assumed to be independent of the angle (at least over some range of angles). A motivation for 
this assumption is that a different distribution of the ejected mass affects only evolution of the jet 
at very early times, but at later times the jet expands based on the Blandford-McKee's (1976) 
self-similar solution during the relativistic phase and the Sedov- Taylor self-similar solution during 
the non-relativistic phase. Thus, the initial Lorentz factor at angle 6 can be written as 

7o() i (*)(*r+i «*<»<«.. (2) 

The evolution of the Lorentz factor 7(0, r) can be calculated from the following equation 

M(r) 1 2 + M 0l = M(r)+M + e(6)/c 2 , (3) 

where M(r) is the mass of the swept-up medium per unit solid angle. Equation (3) expresses 
the conservation of energy and it applies to only the case without heating of the ejected material 
by the reverse shock. In this paper we neglect sideways expansion of the jet. Evidence for this 
consideration is that an observation of GRB 990123 on 7 February 2000 by HST (Fruchter et 
al. 2000), combined with the previous observations, indicates that the steepening of the R-band 
afterglow light curve of this burst during a period of ~ 378 days is still roughly consistent with 
the edge effect of an ultra-relativistic jet of a fixed opening angle expanding in a homogeneous 
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medium proposed by Meszaros & Rees (1999). Further evidence is that, as argued by Jaunsen et 
al. (2000), the rapid initial decline, the sharp break in the optical light curve, and the spectral 
properties of the GRB 980519 afterglow are best interpreted as being due to an ultra-relativistic 
jet of a fixed opening angle expanding in a wind medium with n oc r~ 2 . 

The external medium density is assumed to be a power-law function of radius: n(r) = Ar~ s , 
where A = x lcm~ 3 for an s = (homogeneous) medium (e.g., ISM), A = 3 x lO 35 ^* cm -1 
for an s = 2 wind possibly from a Wolf-Rayet star as the GRB progenitor (Chevalier & Li 1999, 
2000), and A = lO 30 ^* cm -3 / 2 for an s = 3/2 wind possibly from a typical red supergiant star 
(Fransson, Lundquist & Chevalier 1996). For convenience, x is defined as the radius r scaled to 



(3 
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if s = 2 
if s = 3/2, 



(4) 



where £53 = £(#)/10 53 ergs and 70,2 = 7o(#)/10 2 . Furthermore, we have 



= m p n(r)r\ 
3 — s 



(5) 



where m p is the proton mass. 

The solution of equation (3) is 



7(0, r) 



1 

2M 



^M 2 + 4M(M + M + e/c 2 ) - M 



(6) 



The dimensionless radius x and the time t' in the frame comoving with the jet evolve with the 
observer time t based on 



dx c(3 
~dt = r a {l-[ifiy 

M_ 1 
~dt 



(7) 



7(1-/3/")' 

where (5 = yjl — I/7 2 , n = cos 6*, and the line of sight has been taken to be the jet axis. The 
solution of equations (7) and (8) combined with (6) gives the dynamics of a sub-jet at angle 
within the jet. 



(8) 



2.2. Synchrotron Radiation 

We consider synchrotron radiation of the electrons accelerated by the shock. To calculate 
the spectrum and light curve, one needs to determine three crucial frequencies: the synchrotron 
self-absorption frequency (v' a ), the typical synchrotron frequency (v' m ) and the cooling 
frequency (v' c ). We assume a power law distribution of the electrons accelerated by the shock: 
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dn e /d^ e oc % p for 7 e > j m , where p is the index of the electron energy distribution and 
7m = [(p — 2)/(p — l)](m p /m e )e e j is the minimum Lorentz factor. If €b is assumed to be the ratio 
of the magnetic to thermal energy densities of the shocked medium, the magnetic field strength in 
the shocked medium can be approximated by B' = [327re#7(7 — l)nmpC 2 ] 1 / 2 . 

Under these assumptions, the typical synchrotron frequency can be written: 
v' m = 7 m eS' /(2nm e c). The cooling frequency, v' c = ^ 2 eB' /(27rm e c), is clearly determined 
by the Lorentz factor j c , at which an electron is cooling on the dynamical expansion 
timescale (t', measured in the comoving frame). According to Sari, Piran & Narayan (1998), 
7 C = 67rm e c/(<JTB' 2 t f ) where cjt is the Thomson scattering cross-section. Scaled with the typical 
values of the involved parameters, these two frequencies, measured in the comoving frame, are 
derived, 
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vOOlV c 53 ^* 7o,2' 

where t' is in units of 1 s. The synchrotron self-absorption frequency is usually smaller than v' m and 
v' c a few hours after the burst. In this case, the optical depth due to synchrotron self-absorption 
at v' a can be approximated by 

r / /\-5/3 

, 5 enr / v ' \ 

Tab{v' a ) ~ ir- 77 = !. ( n ) 



sfi'75 



where 7 P = min(7 m ,7 c ) and zZ = min(f^ n ,f c ) (see, e.g., Panaitescu & Kumar 2000). At early 
times of the afterglow, all the electrons accelerated behind the shock cool in a timescale smaller 
than the dynamical expansion timescale of the jet at angle 6, i.e., they are in the fast cooling 
regime, implying u' c < v' m , in which case we find the self-absorption frequency, 
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At later times, > f c , in which case the self-absorption frequency is given by 
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The synchrotron peak specific luminosity of a ring (6 
given by 



9 + d9) in the comoving frame is 



dL' m = I 



6.1 x 10 29 ergs s _1 Hz 



-1 I e B U/2 



vO.Ol. 
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In the case of v' c < v' m , the spectrum is (Sari et al. 1998) 



dL' u , 
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If ^ > f^, the spectrum becomes (Sari et al. 1998) 
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(14) 



(15) 
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The observed total flux density at observed frequency v is given by 



dL' 



^7(1 -/3m) 



4vr J D 2 7 3( 1 _ / 3 /U )3' 



(17) 



where Dl is the luminosity distance to the source. For a flat universe with 

Hq = 65kms~ 1 Mpc -1 , the luminosity distance to the source D L = 2c/H (l + z- ^/T+~z) = 

1.8 x 10 28 cm[(l + z)^] 1 / 2 [{VTTz~ - 1)/(V2 - 1)]. 



3. Numerical Results 

Integrating equations (7), (8), and (17) numerically, we can easily obtain an afterglow light 
curve. Figures 1-3 exhibit different afterglow light curves at R-band [y = 4.4 x 10 14 Hz) for 
different k in the homogeneous ISM (n* = 1), in the wind medium with s = 2 (A* = 1), and in 
the wind medium with s = 3/2 (A* = 0.01), respectively. We choose the remaining parameters: 
e = 10 53 ergs sr" 1 , M = 5 x lO^A^sr -1 , 6 = M c 2 /(e + M c 2 ) » 10~ 2 rad, 6 m = 0.3 rad, 
p = 2.5, e e = 0.1, €b = 0.01, and z = 1. In each figure, the lines from top to bottom correspond to 
k = 0, 0.5, 1, 2 and 3 respectively. From these figures, we can see the following features: 

1. At very early times, the afterglow from a jet expanding in the ISM brightens rapidly but the 
afterglow flux for the wind case is approximately a constant, no matter whether the jet is 
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isotropic or anisotropic. This may provide a way of distinguishing between the environments 
as well as GRB progenitors in the near future because rapid, accurate locations of HETE-2 
will be able to allow follow-up observations on very early-time X-ray and optical afterglows. 

2. For an isotropic jet (k = 0) expanding in the ISM, there is a sharp break in the light curve 
at later times. The time at which this break occurs is close to the analytical result, 

f jet « 2 x 10 6 (e 5 3K) 1/3 (# m /0.3) 8 / 3 s, (18) 

at which the Lorentz factor of the jet is equal to the inverse of its half opening angle 
(Meszaros & Rees 1999). But for an isotropic jet expanding in the wind, the break of the 
light curve is weaker and smoother. 

3. In each type of medium, the break of the light curve becomes weaker and smoother as k 
increases. When k = 2 or 3, the break seems to disappear but the afterglow decays rapidly. 
This result shows that the emission from expanding, highly anisotropic jets may provide an 
explanation for some rapidly fading afteglows whose light curves have no break (e.g., GRB 
991208; see the next section). We note that in the preliminary analysis of Meszaros et al. 
(1998), an afterglow decays as one single power-law for any k. The reason for this difference 
will be discussed in the final section. 

4. At very late times, the slopes of the light curves are approximately equal. We also find that 
the slopes become steeper with the electron distribution index p. 

In Figures 4-6 we present the observed radiation spectra computed for a sequence of time for 
a jet with k = 3 in three kinds of medium for the same parameters as in Figures 1-3. It can be 
seen that even for such an anisotropic jet considered here, the observed radiation spectra can still 
be divided into four parts, which can be well described from low-energy to high-energy bands by 
power-law functions F u oc v& , with (3 = 2, 1/3, -(p- l)/2 = -0.75 and -p/2 = -1.25 at t > 10 4 s, 
respectively. Two adjacent spectrum portions are still joined very sharply. These results may be 
due to that the observed radiation spectra are mainly contributed by the < 6$ part of the jet. 

4. The Afterglow of GRB 991208 

GRB 991208 was a long burst with a duration of ~ 60 s and a fluence of ~ 10 -4 ergcm~ 2 
(> 25 keV) (Hurley et al. 1999). Its redshift was measured as z = 0.7055 ± 0.0005 (Djorgovski et 
al. 1999), and thus its isotropic energy in 7-rays is estimated: Ei so ~ 1.3 x 10 53 ergs. Because its 
afterglows at optical, millimeter and radio wavelengths were detected and unusually bright (Sagar 
et al. 2000a; Castro-Tirado et al. 2000; Galama et al. 2000), a detailed study of this burst may 
provide new clues regarding the origin of the GRB phenomenon. 

Sagar et al. (2000a) gave the observed spectral and temporal indexes of the optical afterglow 
/3 G b = —0.75 ±0.03 and a Q b = —2.2 ±0.1 while Galama et al. (2000) presented well-sampled spectra 
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and light curves between radio and millimeter wavelengths for a two-week period, and obtained, 
for the first time, the evolution of the synchrotron self-absorption frequency v a oc 

r 0.15±0.1 2> the 

peak frequency u m oc t _1 - 7±0 - 4 , and the peak flux density F m oc t- - 47 ^ - 11 . The existence of one 
single power law decay for the R-band afterglow implies that the afterglow might be produced by 
one of the following models: (i) an isotropic jet with sideways expansion, (ii) a relativistic fireball 
expanding in the ISM, and (Hi) a relativistic fireball expanding in the wind. Galama et al. (2000) 
argued that model (i) can explain the observations, but models (ii) and (Hi) are ruled out from 
a oh and @ oh . 

If p = 2.5, inferred from /? G b = — (p — l)/2 = —0.75 ± 0.03, Figures 7 and 8 provide two good 
fits to the observed R-band afterglow ligh curve based on our jet model with k = 3 in cases of 
media (ISM and wind with s = 2), respectively. This shows that GRB 991208 might arise from 
a highly anisotropic jet. However, we note that these two cases may not provide further fits to 
the observed radiation spectra because the evolution of the theoretical synchrotron self-absorption 
frequency, peak frequency and peak flux density is inconsistent with the observations by Galama 
et al. (2000). 

Now we assume a generic wind case: n oc r~ s . In this case, we have derived the synchrotron 
self-absorption frequency v a oc t~ 3s /[ 5 ( 4 ~ s )] if p a < i> m , and the peak frequency v m oc t~ 3 / 2 , and the 
peak flux density F m oc i~ s /( 8 ~ 2s ) in the ultra-relativistic phase (Dai & Lu 1998). If s = 3/2, we 
easily find v a oc i~ 0,36 , u m oc t~ 1,5 , and F m oc t~ ' 3 . These scaling laws are in approximate accord 
with those obtained from the observations. This preliminary result prompts us to re-consider a 
highly anisotropic jet expanding in an s = 3/2 wind to fit the afterglow data of GRB 991208. In 
Figure 9 we first show a satisfactory fit to the observed R-band light curve. It is interesting to note 
that such a wind has been inferred in the circumstellar medium of SN 1993J (whose progenitor is 
a red supergiant) by Fransson et al. (1996), and is possibly caused by a variation of the mass- loss 
rate from the progenitor or by a non-spherical geometry. In Figure 10 we further present a fit to 
the radio-to-optical-band afterglow spectrum observed on 1999 December 15.5 UT for the same 
parameters as in Figure 9. Therefore, we can conclude that GRB 991208 might come from a 
highly anisotropic jet expanding in the wind environment (s = 3/2) from a red supergiant. 

5. Discussion and Conclusions 

The model presented in section 2 clearly includes the following assumptions. First, the initial 
opening angle of a highly anisotropic jet (e.g., k = 3), 9 m , was taken to be 0.3, and the observer's 
angle (# bs) between the axis line of the jet and the line of sight to be zero. Because dL m oc e(9) 
in the ultra-relativistic limit, the emission from the jet mainly arises from the shock- accelerated 
electrons moving along the axis line. Thus, the flux density of the emission is weakly dependent 
of 9 m , but strongly depends on 6> b s - For example, the flux density at 6> b s = 0.3 is about four 
orders of magnitude smaller than that of the emission at # b s = 0. GRB 980425 was likely such 
an off-axis burst surrounded by an s = 2 wind medium (Nakamura 1999) because its gamma-ray 
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luminosity of GRB 980425 was also approximately four orders of magnitude smaller than that of 
GRB 970228/980326 and its X-ray afterglow slowly declined. Second, we considered jets of a fixed 
opening angle and neglected the effect of sideways expansion based on two observational facts. 
As shown analytically by Rhoads (1999) and Sari, Piran & Halpern (1999), this effect can lead 
to a sharp break in an afterglow light curve. However, numerical studies by many authors (e.g, 
Panaitescu & Meszaros 1999; Moderski, Sikora & Bulik 2000; Huang et al. 2000a, b; Kumar & 
Panaitescu 2000; Gou et al. 2000; Wei & Lu 2000) show that the actual break, when two effects 
such as the equal-time surface and detailed dynamics of the jet are considered, is much weaker and 
smoother than that predicted analytically. Finally, we used the adiabatic solution for kinematics 
of the jet. This argument is based on the radiative efficiency of the jet, / = e e [*syn/(*syn + *cxp)]> 
defined by Dai et al. (1999). At very early times, the cooling timescale due to synchrotron 
radiation, t' syn , may be much less than the expansion timescale, tg xp , and / f» 0.1(e e /0.1); but at 
late times, t' syn » t' cxp and thus / « 0.1(e e /0.1)(ig Xp /ig yn ) <C 1. Therefore, it can be seen that the 
energy losses due to synchrotron radiation are insignificant during the whole evolution stage of the 
jet. 

Our numerical results show that for an isotropic jet (k = 0) expanding in the ISM, there is 
a sharp break in the light curve at late times, but for an isotropic jet expanding in the wind, the 
break of the light curve is weaker and smoother. This result is easily understood in terms of the 
analytical model of Dai & Lu (1998) and Chevalier & Li (1999, 2000): at early times of the jet 
evolution, 7 > in which stage the jet is spherical-like and thus the decay index of the afterglow 
a.\ = — [s/(8 — 2s) + 3(p — l)/4] in the slow cooling regime; but at later times, 7 < O^ 1 , the decay 
index becomes 02 = — [(6 — s)/(8 — 2s) + 3(p — l)/4] due to the edge effect. These indexes imply 
the steepening of the light curve Aa = a\ — Q2 = (3 — s)/(4 — s) = 0.5 for s = 2, 0.6 for s = 3/2, 
or 0.75 for s = 0. Therefore, the break in the light curve of the afterglow from a jet expanding in 
a wind is weaker and smoother than that for a homogeneous medium case. Our results also show 
that for an anisotropic jet, one break seems to appear in the light curve at small k but becomes 
weaker and smoother as k increases. This conclusion is different from that presented by Meszaros 
et al. (1998), who found that an afterglow decays as one single power-law for any k. The reason 
for this difference is that we took into account the light travel effects related to different sub-jets 
within the jet but Meszaros et al. (1998) didn't. It is easily understood that an anisotropic jet 
with a small k can be treated as a quasi-isotropic jet with a half opening angle 8 m , and thus the 
temporal decay of the flux density of radiation from such a jet will begin to expedite at t = tj e t, 
which is given by equation (18) for the ISM case. However, for a highly anisotropic jet (e.g., k = 2 
or 3), the contribution of radiation from 8 > 8q is in fact much smaller than that from 8 < 8q due 
to the fact that the energy per unit solid angle (dE/dQ) within this jet decreases rapidly with 
increasing 8, and thus, the light curve of the resulting radiation should begin to steepen at a very 
early time, t « 230(e53/n*) 1 / 3 (^o/10 _2 ) 8//3 s, which is also estimated from equation (18) in the 
ISM case but for 8 m being replaced by 8q. Therefore, one cannot see a broken light curve at late 
times for such a highly anisotropic jet. 
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In summary, the energy per unit solid angle within a realistic jet may be a power-law function 
of the angle (oc 9~ k ). Such anisotropic jets are expected in the collapsar model of MacFadyen et 
al. (2000). We numerically calculated light curves and spectra of the emission from anisotropic 
jets expanding either in the interstellar medium (ISM) or in the wind medium. We also took into 
account two kinds of wind: one (n oc r~ 3 / 2 ) possibly from a typical red supergiant and another 
(n oc r -2 ) possibly from a Wolf-Rayet star. Two of the main results of this work are that (i) at very 
early times, the afterglow from a jet expanding in the ISM brightens rapidly but the afterglow flux 
for the wind case is approximately a constant, no matter whether the jet is isotropic or anisotropic. 
Based on this conclusion, future observations led by HETE-2 on very early-time optical afterglows 
may be able to distinguish between the environments as well as GRB progenitors, (it) In each 
type of medium, one break appears in the late-time afterglow light curve for small k, but becomes 
weaker and smoother as k increases. When k > 2, the break seems to disappear but the afterglow 
decays rapidly. Therefore, one expects that the emission from expanding, highly anisotropic jets 
provides an explanation for some rapidly fading afteglows whose light curves have no break. We 
also presented good fits to the optical afterglow flux data of GRB 991208. Finally, to interpret the 
observed radio-to-optical-band afterglow data of this burst, we argued that it might arise from a 
highly anisotropic jet expanding in the wind (n oc r~ 3 / 2 ) from a red supergiant. 
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Fig. 1. — Light curves of R-band afterglows from anisotropic jets expanding in the ISM (s = 0). The 
lines from top to bottom correspond to k = 0, 0.5, 1, 2 and 3 respectively. The other parameters 
for the model are seen in the text. 
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Fig. 2. — The same as in Figure 1 but for the s = 2 wind. 
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Fig. 3. — The same as in Figure 1 but for the s = 3/2 wind. 
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Fig. 4. — Evolution of the observed radiation spectra for a k = 3 jet expanding in the ISM (s = 0). 
The numbers mark the observed times in units of second. The other parameters for the model are 
the same as in Figure 1. 



-17- 




Fig. 5 — 



The same as in Figure 4 but for the s = 2 wind. 
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Fig. 6. — The same as in Figure 4 but for the s = 3/2 wind. 
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Fig. 7. — Comparison between the observed and theoretically calculated light curves for the ri- 
band afterglow of GRB 991208. The data are taken from Sagar et al. (2000a) and Castro-Tirado 
et al. (2000), and the model light curve is calculated for a k = 3 jet expanding in the ISM (s = 0) 



when an observer is located on the jet axis. The model parameters are chosen: eq = 10 53 ergs sr 



M = 5 x lO- 4 M sr-\ 9 = 10" 
The insert shows a clearer fit. 



rad, 9 m = 0.3 rad, p = 2.5, = 1.0, e e = 0.25 and e B = 0.01. 
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Fig. 8. — The same as in Figure 7 but for s = 2 (wind), A* = 



1.0 and e e = 0.21. 
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Fig. 9. — The same as in Figure 7 but for the s = 3/2 wind case and for the following model 
parameters: e = 1-5 x 10 53 ergs sr^ 1 , 6 = 0.01 rad, A* = 8 x 10~ 3 , e e = 0.16, e# = 10~ 3 , p = 2.5, 
k = 3, and m = 0.3 rad. 
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Fig. 10. — Comparison between the observed and theoretically calculated spectra of the GRB 
991208 afterglow on 1999 December 15.5 UT. The data are taken from Galama et al. (2000), Sagar 
et al. (2000a) and Castro-Tirado et al. (2000). The model and its parameters are the same as in 
Figure 9. 



